Strong selection signatures for Aleutian disease tolerance acting on novel candidate genes linked to immune and cellular responses in American mink (Neogale vison).

Aleutian disease (AD) is a multi-systemic infectious disease in American mink (Neogale vison) caused by Aleutian mink disease virus (AMDV). This study aimed to identify candidate regions and genes underlying selection for response against AMDV using whole-genome sequence (WGS) data. Three case-control selection signatures studies were conducted between animals (N = 85) producing high versus low antibody levels against AMDV, grouped by counter immunoelectrophoresis (CIEP) test and two enzyme-linked immunosorbent assays (ELISA). Within each study, selection signals were detected using fixation index (FST) and nucleotide diversity (θπ ratios), and validated by cross-population extended haplotype homozygosity (XP-EHH) test. Within- and between-studies overlapping results were then evaluated. Within-studies overlapping results indicated novel candidate genes related to immune and cellular responses (e.g., TAP2, RAB32), respiratory system function (e.g., SPEF2, R3HCC1L), and reproduction system function (e.g., HSF2, CFAP206) in other species. Between-studies overlapping results identified three large segments under strong selection pressure, including two on chromosome 1 (chr1:88,770-98,281 kb and chr1:114,133-120,473) and one on chromosome 6 (chr6:37,953-44,279 kb). Within regions with strong signals, we found novel candidate genes involved in immune and cellular responses (e.g., homologous MHC class II genes, ITPR3, VPS52) in other species. Our study brings new insights into candidate regions and genes controlling AD response.

FSTest: an efficient tool for cross-population fixation index estimation on variant call format files.

Fixation index (Fst) statistics provide critical insights into evolutionary processes affecting the structure of genetic variation within and among populations. Fst statistics have been widely applied in population and evolutionary genetics to identify genomic regions targeted by selection pressures. The FSTest 1.3 software was developed to estimate four Fst statistics of Hudson, Weir and Cockerham, Nei, and Wright using high-throughput genotyping or sequencing data. Here, we introduced FSTest 1.3 and compared its performance with two widely used software VCFtools 0.1.16 and PLINK 2.0. Chromosome 1 of 1000 Genomes Phase III variant data belonging to South Asian (n = 211) and African (n = 274) populations were included as an example case in this study. Different Fst estimates were calculated for each single-nucleotide polymorphism (SNP) in a pairwise comparison of South Asian against African populations, and the results of FSTest 1.3 were confirmed by VCFtools 0.1.16 and PLINK 2.0. Two different sliding window approaches, one based on a fixed number of SNPs and another based on a fixed number of base pair (bp) were conducted using FSTest 1.3 and VCFtools 0.1.16. Our results showed that regions with low coverage genotypic data could lead to an overestimation of Fst in sliding window analysis using a fixed number of bp. FSTest 1.3 could mitigate this challenge by estimating the average of consecutive SNPs along the chromosome. FSTest 1.3 allows direct analysis of VCF files with a small amount of code and can calculate Fst estimates on a desktop computer for more than a million SNPs in a few minutes. FSTest 1.3 is freely available at https://github.com/similab/FSTest.

Population genetic analysis and scans for adaptation and contemporary selection footprints provide genomic insight into aus, indica and japonica rice cultivars diversification.

Following domestication, rice cultivars have been spread worldwide to different climates and have experienced selection pressures to improve desirable traits. This has resulted in diverse cultivars that display variations in phenotypic traits, such as stress tolerance, grain size, and yield. To better understand the genomic composition arising from cultivar's development and local adaptation, high-density genotypes (containing 286,183 single-nucleotide polymorphisms after the quality control) of 1284 rice cultivars of aus, indica, and temperate and tropical japonica were scanned for diversifying signatures by applying a pairwise comparison of fixation index (Fst) test. Each cultivar's population was investigated for contemporary selection using the integrated haplotype score test. Signatures of diversifying selection among the pairwise comparisons were found in genomic regions mainly involved in response to stress (pathogens, drought, heat, cold) and development and morphology of various structures, such as root, pollen, spikelet, and grain. The most significant diversification signal between indica and japonica cultivars was detected at the location of ROX2 gene. Aus with indica comparison detected the most divergent signal at important candidate genes of OsEXPA8 and OsEXPA9, whereas temperate with tropical japonica comparison resulted in two well-known candidate genes OsHCT4 and OsGpx4. Recent selection analysis detected different patterns of contemporary selection in genomic regions related to rice breeding standard criteria such as stress tolerance, seed germination, starch content, and flowering time. Our findings highlight the underlying molecular basis of adaptive divergence and propose that modern rice breeding may provide additional diversification among rice cultivars.

Epidemiology, pathogenesis, and diagnosis of Aleutian disease caused by Aleutian mink disease virus: A literature review with a perspective of genomic breeding for disease control in American mink (Neogale vison).

Aleutian disease (AD) is a multi-systemic infectious disease in American mink (Neogale vison) caused by the Aleutian mink disease virus (AMDV). Commonly referred to as mink plasmacytosis, AD is an economically significant disease in mink-breeding countries. Aleutian disease mainly induces weight loss, lower fertility, and dropped pelt quality in adults and can result in acute interstitial pneumonia with high mortality rates in kits. In this review, we employed the scientific literature on AD over the last 70 years to discuss the historical and contemporary status of AD outbreaks and seroprevalence in mink farming countries. We also explained different forms of AD and the differences between the pathogenicity of the virus in kits and adults. The application of the available AD serological tests in AD control strategies was argued. We explained how selection programs could help AD control and proposed different approaches to selecting animals for building AD-tolerant herds. The advantages of genomic selection for AD tolerance over traditional breeding strategies were discussed in detail. We also explained how genomic selection could help AD control by selecting tolerant animals for the next generation based on genome-wide single nucleotide polymorphisms (SNP) data and the challenges of implementing genomic selection for AD tolerance in the mink industry. This review collected the information required for designing successful breeding programs for AD tolerance. Examples of the application of information are presented, and data gaps are highlighted. We showed that AD tolerance is necessary to be among the traits that animals are selected for in the mink industry.

Expanding the application of haplotype-based genomic predictions to the wild: A case of antibody response against Teladorsagia circumcincta in Soay sheep.

BACKGROUND: Genomic prediction of breeding values (GP) has been adopted in evolutionary genomic studies to uncover microevolutionary processes of wild populations or improve captive breeding strategies. While recent evolutionary studies applied GP with individual single nucleotide polymorphism (SNP), haplotype-based GP could outperform individual SNP predictions through better capturing the linkage disequilibrium (LD) between the SNP and quantitative trait loci (QTL). This study aimed to evaluate the accuracy and bias of haplotype-based GP of immunoglobulin (Ig) A (IgA), IgE, and IgG against Teladorsagia circumcincta in lambs of an unmanaged sheep population (Soay breed) based on Genomic Best Linear Unbiased Prediction (GBLUP) and five Bayesian [BayesA, BayesB, BayesCπ, Bayesian Lasso (BayesL), and BayesR] methods. RESULTS: The accuracy and bias of GPs using SNP, haplotypic pseudo-SNP from blocks with different LD thresholds (0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.00), or the combinations of pseudo-SNPs and non-LD clustered SNPs were obtained. Across methods and marker sets, higher ranges of genomic estimated breeding values (GEBV) accuracies were observed for IgA (0.20 to 0.49), followed by IgE (0.08 to 0.20) and IgG (0.05 to 0.14). Considering the methods evaluated, up to 8% gains in GP accuracy of IgG were achieved using pseudo-SNPs compared to SNPs. Up to 3% gain in GP accuracy for IgA was also obtained using the combinations of the pseudo-SNPs with non-clustered SNPs in comparison to fitting individual SNP. No improvement in GP accuracy of IgE was observed using haplotypic pseudo-SNPs or their combination with non-clustered SNPs compared to individual SNP. Bayesian methods outperformed GBLUP for all traits. Most scenarios yielded lower accuracies for all traits with an increased LD threshold. GP models using haplotypic pseudo-SNPs predicted less-biased GEBVs mainly for IgG. For this trait, lower bias was observed with higher LD thresholds, whereas no distinct trend was observed for other traits with changes in LD. CONCLUSIONS: Haplotype information improves GP performance of anti-helminthic antibody traits of IgA and IgG compared to fitting individual SNP. The observed gains in the predictive performances indicate that haplotype-based methods could benefit GP of some traits in wild animal populations.

Genome-wide assessment and mapping of inbreeding depression identifies candidate genes associated with semen traits in Holstein bulls.

BACKGROUND: The reduction in phenotypic performance of a population due to mating between close relatives is called inbreeding depression. The genetic background of inbreeding depression for semen traits is poorly understood. Thus, the objectives were to estimate the effect of inbreeding and to identify genomic regions underlying inbreeding depression of semen traits including ejaculate volume (EV), sperm concentration (SC), and sperm motility (SM). The dataset comprised ~ 330 K semen records from ~ 1.5 K Holstein bulls genotyped with 50 K single nucleotide polymorphism (SNP) BeadChip. Genomic inbreeding coefficients were estimated using runs of homozygosity (i.e., FROH > 1 Mb) and excess of SNP homozygosity (FSNP). The effect of inbreeding was estimated by regressing phenotypes of semen traits on inbreeding coefficients. Associated variants with inbreeding depression were also detected by regressing phenotypes on ROH state of the variants. RESULTS: Significant inbreeding depression was observed for SC and SM (p < 0.01). A 1% increase in FROH reduced SM and SC by 0.28% and 0.42% of the population mean, respectively. By splitting FROH into different lengths, we found significant reduction in SC and SM due to longer ROH, which is indicative of more recent inbreeding. A genome-wide association study revealed two signals positioned on BTA 8 associated with inbreeding depression of SC (p < 0.00001; FDR < 0.02). Three candidate genes of GALNTL6, HMGB2, and ADAM29, located in these regions, have established and conserved connections with reproduction and/or male fertility. Moreover, six genomic regions on BTA 3, 9, 21 and 28 were associated with SM (p < 0.0001; FDR < 0.08). These genomic regions contained genes including PRMT6, SCAPER, EDC3, and LIN28B with established connections to spermatogenesis or fertility. CONCLUSIONS: Inbreeding depression adversely affects SC and SM, with evidence that longer ROH, or more recent inbreeding, being especially detrimental. There are genomic regions associated with semen traits that seems to be especially sensitive to homozygosity, and evidence to support some from other studies. Breeding companies may wish to consider avoiding homozygosity in these regions for potential artificial insemination sires.

Population structure and genomic footprints of selection in five major Iranian horse breeds.

The genetic structure and characteristics of Iranian native breeds are yet to be comprehensibly investigated and studied. Therefore, we employed genomic information of 364 Iranian native horses representing the Asil (n = 109), Caspian (n = 40), Dareshuri (n = 44), Kurdish (n = 95), and Turkoman (n = 76) breeds to reveal the genetic structure and characteristics. For these and 19 other horse breeds, principal component analysis, Bayesian model-based, Neighbor-Net, and bootstrap-based TreeMix approaches were applied to investigate and compare their genetic structure. Additionally, three haplotype-based methods including haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length were applied to trace genomic footprints of selection of Asil, Caspian, Dareshuri, Kurdish, and Turkoman groups. Then, the Mahalanobis distance based on the negative-log10 rank-based P-values was estimated based on the haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length values. Asil, Caspian, Dareshuri, Kurdish, and Turkoman can be categorized into five different genetic clusters. Based on the top 1% of Mahalanobis distance based on the negative-log10 rank-based P-values of SNPs, we identified 24 SNPs formerly reported to be associated with different traits and >100 genes undergoing selection pressures in Asil, Caspian, Dareshuri, Kurdish, and Turkoman. The detected QTL undergoing selection pressures were associated with withers height, equine metabolic syndrome, overall body size, insect bite hypersensitivity, guttural pouch tympany, white markings, Rhodococcus equi infection, jumping test score, alternate gaits, and body weight traits. Our findings will aid to have a better perspective of the genetic characteristics and population structure of Asil, Caspian, Dareshuri, Kurdish, and Turkoman horses as Iranian native horse breeds.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle.

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The fdM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

Detection of Common Copy Number of Variants Underlying Selection Pressure in Middle Eastern Horse Breeds Using Whole-Genome Sequence Data.

Dareshouri, Arabian, and Akhal-Teke are 3 Middle Eastern horse breeds that have been selected for endurance and adaptation to harsh climates. Deciphering the genetic characteristics of these horses by tracing selection footprints and copy number of variations will be helpful in improving our understanding of equine breeds' development and adaptation. For this purpose, we sequenced the whole genome of 4 Dareshouri horses using Illumina Hiseq panels and compared them with publicly available whole-genome sequences of Arabian (n = 3) and Akhal-Teke (n = 3) horses. Three tests of FLK, hapFLK, and pooled heterozygosity were applied using a sliding window (window size = 100 kb, step size = 50 kb) approach to detect putative selection signals. Copy number variation analysis was applied to investigate copy number of variants (CNVs), and the results were used to suggest selection signatures involving CNVs. Whole-genome sequencing demonstrated 8 837 950 single-nucleotide polymorphisms (SNPs) in autosomal chromosomes. We suggested 58 genes and 3 quantitative trait loci, including some related to horse gait, insect bite hypersensitivity, and withers height, based on selective signals detected by adjusted P-value of Mahalanobis distance based on the rank-based P-values (Md-rank-P) method. We proposed 12 genomic regions under selection pressure involving CNVs that were previously reported to be associated with metabolism energy (SLC5A8), champagne dilution in horses (SLC36A1), and synthesis of polyunsaturated fatty acids (FAT2). Only 10 Middle Eastern horses were tested in this study; therefore, the conclusions are speculative. Our findings are useful to better understanding the evolution and adaptation of Middle Eastern horse breeds.

Weighted Single-Step GWAS for Body Mass Index and Scans for Recent Signatures of Selection in Yorkshire Pigs.

Controlling extra fat deposition is economically favorable in modern swine industry. Understanding the genetic architecture of fat deposition traits such as body mass index (BMI) can help in improving genomic selection for such traits. We utilized a weighted single-step genome-wide association study (WssGWAS) to detect genetic regions and candidate genes associated with BMI in a Yorkshire pig population. Three extended haplotype homozygosity (EHH)-related statistics were also incorporated within a de-correlated composite of multiple signals (DCMS) framework to detect recent selection signatures signals. Overall, the full pedigree consisted of 7016 pigs, of which 5561 had BMI records and 598 pigs were genotyped with an 80 K single nucleotide polymorphism (SNP) array. Results showed that the most significant windows (top 15) explained 9.35% of BMI genetic variance. Several genes were detected in regions previously associated with pig fat deposition traits and treated as potential candidate genes for BMI in Yorkshire pigs: FTMT, SRFBP1, KHDRBS3, FOXG1, SOD3, LRRC32, TSKU, ACER3, B3GNT6, CCDC201, ADCY1, RAMP3, TBRG4, CCM2. Signature of selection analysis revealed multiple candidate genes previously associated with various economic traits. However, BMI genetic variance explained by regions under selection pressure was minimal (1.31%). In conclusion, candidate genes associated with Yorkshire pigs' BMI trait were identified using WssGWAS. Gene enrichment analysis indicated that the identified candidate genes were enriched in the insulin secretion pathway. We anticipate that these results further advance our understanding of the genetic architecture of BMI in Yorkshire pigs and provide information for genomic selection for fat deposition in this breed.

Intestinal microsporidia infection among cat owners and non-pet owners in Iran: a case-control study.

Microsporidia is a group of spore-forming microorganisms with zoonotic potential. This study aimed to compare intestinal microsporidia infections in cat owners and non-pet owners. In total, 210 fecal samples were collected from indoor cats, cat owners, and non-pet owners. DNA extraction was performed and the small subunit ribosomal RNA (SSU rRNA) gene was amplified. To characterize the genotypes, the internal transcribed spacer (ITS) fragment was amplified and sequenced. The phylogenetic trees were drawn to evaluate the relationship among Enterocytozoon bieneusi isolates. Two (2.9%) and one (1.4%) fecal samples from cat owners and one (1.4%) and two (2.9%) fecal samples from non-pet owners were positive for E. bieneusi and Encephalitozoon intestinalis, respectively. E. bieneusi was detected in two cat samples (2.9%). Same infection was not seen between infected cats and their owners. There was no significant difference between the prevalence rate of microsporidia among the cat owners and non-pet owners. Indeed, the genotypes L and type IV were seen in cats, while the genotype D was only detected in human. In this study, E. bieneusi and E. intestinalis were more prevalent among the cat owners and non-pet owners, respectively. Indeed, the higher prevalence of E. bieneusi in cats and their owners might be resulted from the worldwide distribution of this species.

Antiviral responses against chicken respiratory infections: Focus on avian influenza virus and infectious bronchitis virus.

Some of the respiratory viral infections in chickens pose a significant threat to the poultry industry and public health. In response to viral infections, host innate responses provide the first line of defense against viruses, which often act even before the establishment of the infection. Host cells sense the presence of viral components through germinal encoded pattern recognition receptors (PRRs). The engagement of PRRs with pathogen-associated molecular patterns leads to the induction of pro-inflammatory and interferon productions. Induced antiviral responses play a critical role in the outcome of the infections. In order to improve current strategies for control of viral infections or to advance new strategies aimed against viral infections, a deep understanding of host-virus interaction and induction of antiviral responses is required. In this review, we summarized recent progress in understanding innate antiviral responses in chickens with a focus on the avian influenza virus and infectious bronchitis virus.

Association between DLA-DRB1.2 allelic diversity and development of mammary gland tumors in dogs.

BACKGROUND: The major histocompatibility complex (MHC) is the best-characterized genetic region related to resistance/susceptibility to a wide range of infectious and immune-mediated diseases. Evidences suggest that MHC class II genes may play an important role in developing different types of tumors including breast cancer. Canine mammary gland tumors (CMTs) are the most common neoplasms in female dogs. In the current study, the association of canine MHC class II DLA-DRB1.2 genotypes with development of mammary gland tumor profiles in dogs was investigated. DLA-DRB1.2 allelic diversity was determined in 40 dogs (18 CMT cases and 22 controls) using HRM technique and DNA sequencing. Association of the DLA-DRB1.2 genotypes with CMT profiles was expressed as odds ratio (OR). RESULTS: Based on the histopathological typing of tumors, CMT cases were categorized into 4 groups: simple carcinoma, complex carcinoma, carcinoma arising in a benign tumor and special types of carcinoma. A total of eight HRM profiles (A to H) were identified in dogs sampled. The association study revealed a significant correlation between DLA-DRB1.2 genotypes with different CMT profiles. The E genotype was significantly associated with increased risk of carcinoma arising in a benign tumor, and the B genotype represented a positive correlation with complex carcinoma. Significant association was also observed between the heterozygosity of DLA-DRB1.2 genotypes and decreased risk of developing tumor in dogs. CONCLUSIONS: These results provide additional support for the association between DLA-DRB1 genes and development of mammary gland tumors in dogs and could potentially be used for early diagnosis of neoplasia and identifying susceptible dogs.

Molecular identification and phylogenic analysis of Bartonella henselae isolated from Iranian cats based on gltA gene.

One of the most important species of the Bartonella genus is B. henselae that causes a zoonotic infection, cat scratch disease (CSD). The main source of the bacteria is cat and the carrier is Ctenocephalides felis flea. One hundred and forty nail and saliva samples were collected from 70 domestic cats. Positive samples for B. henselae were characterized by polymerase chain reaction (PCR) and sequencing. Sequences of gltA gene were trimmed using BioEdit software and then compared with the sequences of the same gene from B. henselae isolated from cats and humans in GenBank database. Phylogenic tree was constructed using CLC Sequence Viewer software and unweighted pair group method with arithmetic mean (UPGMA) method. Molecular assessments showed that five samples out of 70 nail samples (7.14%) and one sample out of 70 saliva samples (1.42%) were genetically positive for B. henselae. At least an 87.00% similarity was seen between the gene sequences from the current study and the reference sequences from the GenBank database. Phylogenic analysis has shown that strains isolated in this study were grouped in a different haplo group, compared to other strains. Among the Asian countries, the prevalence of the bacteria in Iran was close to that in Japan and Turkey. In conclusion, findings of this study showed the prevalence of B. henselae in Iranian cats which is important due to its public health issues, especially for the immunocompromised pet owners.